Finite-size effects and collective vibrations in the inner crust of neutron stars 
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We study the linear response of the inner crust of neutron stars within the Random Phase Approx- 
imation, employing a Skyrme-type interaction as effective interaction. We adopt the Wigner-Seitz 
approximation, and consider a single unit cell of the Coulomb lattice which constitutes the inner 
crust, with a nucleus at its center, surrounded by a sea of free neutrons. With the use of an appro- 
priate operator, it is possible to analyze in detail the properties of the vibrations of the surface of 
the nucleus and their interaction with the modes of the sea of free neutrons, and to investigate the 
role of shell effects and of resonant states. 
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I. INTRODUCTION 

In going from the outer crust of a neutron star to- 
wards the core of the star, one crosses the so-called inner 
crust. The baryonic density of the inner crust ranges 
from about p — 10 11 gem -3 to approximately po/2, 
where po = 2.8 x 10 14 g cm" 3 (or po = 0.16 fm~ 3 ) is the 
density of nuclear matter at saturation. We shall discuss 
the linear response of this system limiting ourselves to 
densities smaller than po/10. According to the generally 
accepted theoretical description the inner crust consists, 
in this region, of a lattice of spherical nuclei immersed in 
a sea of free neutrons (with a background of electrons) 
P, 0- At these densities, neutrons are expected to be 
superfiuid, and the inner crust of a rotating star should 
be threaded by vortex lines. 

The presence of nuclei affects a number of thermal and 
transport properties of the inner crust, and in some as- 
trophysical scenarios may influence the fast cooling of the 
star [HH}. It can also affect the vortex dynamics, leading 
to vortex pinning, which might be the explanation for the 
phenomenon of glitches, sudden changes observed in the 
rotational period of pulsars 0-Q]- 

The interplay between the degrees of freedom of nuclei 
and those of free neutrons has been considered by differ- 
ent groups at the mean-field level within Hartree-Fock- 
Bogoliubov theory (HFB) 0] , . It is well known, 
however, that many-body effects can significantly mod- 
ify the HFB results. Several studies have found that the 
pairing gap in neutron matter is quenched by spin fluc- 
tuations [ilHl4| . but results obtained in uniform matter 
cannot be directly extrapolated to the actual case of the 
inhomogeneous inner crust. In fact, many-body effects 
in (isolated) atomic nuclei are strongly dominated by the 
coupling between single-particle and collective surface- 
like degrees of freedom. The exchange of surface modes 
between nucleons moving in time-reversal states close to 
the Fermi energy gives rise to an attractive interaction, 
which tends to enhance the pairing gap [l5l flfj| . Only 



exploratory calculations have directly addressed the in- 
duced interaction in the inner crust, finding that the den- 
sity fluctuations associated with the dynamics of the nu- 
clear surface lead to a partial dequenching of the gap and 
to characteristic changes in the spatial dependence of the 
gap H3. 

Within this scenario it is of interest to investigate how 
the surface response changes from atomic nuclei (nega- 
tive value of Fermi energy) to the inner crust, where the 
neutron Fermi energy lies in the continuum. 

In the present work we shall investigate the collective 
response of the nuclear surface by means of a microscopic 
calculation of the linear response of the system based on 
the Random Phase Approximation (RPA) theory filjj . 
As in most previous studies, we shall deal with this inho- 
mogeneous system within the Wigner-Seitz approxima- 
tion (WS), enclosing the system in a spherical unit cell. 
In this way, we shall neglect the band structure associ- 
ated with the Coulomb lattice [l9j . 

Previous studies [20l - [2^ | have shown that the response 
to operators of the type r L YLM or [t l Yl X <j]jm is very 
close to that found in neutron matter, in keeping with 
the fact that the nucleus occupies only a small fraction of 
the volume of the WS cell. However, in the following we 
shall focus on the modes which are specifically associated 
with the presence of the nucleus at the center of the WS 
cell. Therefore, we shall study the response to the oper- 
ator dU /drYLM , where U is the mean-field potential. In 
fact, dU/drY^M is the genuine operator associated with 
collective surface modes, which appears naturally in the 
particle- vibration coupling Hamiltonian [23[ . Essentially 
the same results would be obtained using the density in- 
stead of the mean-field potential. In Appendix B we also 
briefly consider the operator r L Yr,M- 

In this work we shall not consider the effects of super- 
fluidity on the linear response. Pairing correlations can 
affect considerably the low-lying part of the spectrum, 
not only because they smear the Fermi surface leading to 
the partial occupation of single-particle states, but also 
because they can change the isotopic composition of the 



2 



crust [24|. We shall focus on the properties of giant res- 
onances, which are not expected to be sensitive to pair- 
ing. Within this context, we shall use in the calculations 
the favoured number of protons predicted by Negele and 
Vautherin |25[ , predictions which were worked out with- 
out pairing. We have anyway indicated in the text where 
one expects to sec the main effects of pairing correlations. 



II. OUTLINE OF THE CALCULATIONS 

We start from a mean-field Hartree-Fock (HF) calcula- 
tion, making use of the Skyrme-type SLy4 interaction [26j 
and working on a grid with a mesh size of 0.1 fm. Our 
calculation is in line with other studies previously per- 
formed by a number of groups [!, ■ We impose that 
the single-particle wavefunctions vanish at the edges of 
the WS cell. We expect that within the selected density 
range the results should not be sensitive to the particular 
choice of the boundary conditions (cf. also Appendix A) . 
For densities larger than those considered in this paper, 
the WS approximation starts to break down, leading to 
a dependence on the boundary conditions [24j . 

We shall show results for two WS cells. The first 
one has a radius Rws= 33.1 fm and contains 1314 neu- 
trons and 50 protons ( 1364 Sn), corresponding to nucleon 
density p k, 1.5 x 10 13 g cm~ 3 (or 0.05/?o)- Far from 
the nucleus the density approaches the asymptotic den- 
sity poo k, 0.01 fm ~ 3 . The second one has a ra- 
dius Rws = 42.2 fm and contains 458 neutrons and 
40 protons ( 498 Zr), corresponding to the nucleon density 
p w 2.7 x 10 12 g cm" 3 (or 0.01p ) and to p x w 1.5 x 10~ 3 
fm~ 3 . The number of protons and neutrons and the cell 
radius are an input to the HF calculation and have been 
taken from the results of Negele and Vautherin [25[ . They 
found that the favoured proton number is equal to Z = 40 
in the lower-density part of the inner crust, changing to 
Z = 50 for densities larger than p « 3 x 10 12 g cm -3 . 
Other authors have determined the isotopic composition 
of the crust using different energy functionals and dif- 
ferent methods. In particu lar, systematic calculations 
have been performed in ref. [24j including pairing corre- 
lations, neglected by Negele and Vautherin. The results 
of these studies have been summarized in ref. [l[. In 
most cases, one finds that the favoured proton number 
varies between Z m 35 and Z = 50, becoming smaller 
(Z w 25) in the deeper regions of the crust. These uncer- 
tainties are understandable, because typical energy dif- 
ferences between local minima in the energy as a function 
of the proton number at a given density are of the order of 
10 keV/nucleon. In the following, we shall also consider 
the WS cell 506 Sn, calculated with the same parameters 
as 498 Zr, but for the different proton number Z = 50 
instead of Z = 40. 

The corresponding Fermi energy for neutrons and pro- 
tons are equal to Ep = 5.5 MeV and Ep = —35.1 MeV 
for 1364 Sn, and to Ep = 2.1 MeV and Ep = —31.5 MeV 
for 498 Zr. We note that the potential reaches an almost 
constant value far from the nucleus, so that the con- 



tinuum portion of the single-particle spectrum starts at 
about E cont = -3.5 MeV ( 1364 Sn) and at E cont = -0.5 
MeV ( 498 Zr). The number of neutrons in bound orbitals 
is about the same in the 1364 Sn (126 bound neutrons) 
and in the 498 Zr (124 bound neutrons) cases. 

In Fig[T] we show the mean-field potential U(r) in 
the two cases [27|. The presence of the outer neu- 
tron gas leads to a more diffuse potential than that 
found in atomic nuclei. The neutron potentials can be 
parametrized as Woods-Saxon potentials, 



U(r) = U OQ + 



U, 







1 + exp 



r-Ro 



(1) 



with the values U^, = -3.73 MeV, U = -54.2 MeV, 
Ro = 7.06 fm, a = 0.79 fm for 1364 Sn and = -0.42 
MeV, U = -61.4 MeV, Ro = 6.42 fm, a = 0.73 fm for 
498 Zr. We remark that the values of the radii are close to 
those of ordinary nuclei with the same number of bound 
nucleons (1.2 A 1 ^ w 7 fm for A w 200 ). The extension 
of proton and neutron potentials is similar. The radius 
of 498 Zr is close to the value calculated for the drip-line 
nucleus 176 Sn (cf. below Section 111(c)). We can use 
the value of the radius Rq to estimate the typical oscilla- 
tor frequencies associated with the mean-field potentials 
shown in FigfTJ Using the expression Huj « 75V 1 / 3 /i? 2 
(which reduces to the usual estimate hu> « 40^4~ 1 / 3 for 
N = Z and Ro = 1.2A 1 / 3 ), one finds Hu n = 7.7 MeV for 
neutrons (using the number of bound neutrons, N ss 150) 
and hojp = 5.3 MeV for protons in 1364 Sn. Taking into 
account that the effective mass associated with the 
HF potential increases the level spacing, one can estimate 
hw n = 9.6 MeV and hujp = 6.6 MeV, using mk/m = 0.8 
(the typical value of the effective mass associated with 
the SLy4 interaction in the nuclear surface region). We 
note the small value of Hu) p , associated with the low pro- 
ton density. 

We have then calculated the excited states of the sys- 
tem performing a RPA calculation in configuration space, 
based on the HF mean-field. The residual interaction 
was taken as the second derivative of the energy func- 
tional, that is, it was calculated with the same two-body 
Skyrme force which determines the mean-field. The only 
terms that were dropped in the p-h channel were the two- 
body Coulomb and two-body spin-orbit. Their impor- 
tance was recently discussed in refs. [H, [29[ for atomic 
nuclei. While they play a role in a detailed analysis of 
the experimental data, they hardly change the general 
properties of nuclear response. 

The resulting RPA phonons are classified according to 
their total angular momentum J and parity n, while the 
index v = 1,2,... runs over the states having the same 
spin and parity with increasing energy. We have tested 
the numerical scheme considering the case of WS cells 
without protons. In this case, the nuclear potential is flat 
except at the edge of the cell, and the WS cell represents 
a portion of neutron matter. The calculated response can 
then be compared with the analytical results obtained in 
uniform matter. The corresponding results are collected 
in Appendix A. 
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FIG. 1. (left) Mean-field potential calculated in the WS cell 1364 Sn. (right) The same, for the WS cell 498 Zr. The thin horizontal 
lines indicate the neutron and the proton Fermi energies, while in the insets we show the local densities. 



III. RESULTS 

As discussed in the introduction, in the present work 
we are particularly interested in the study of the dynam- 
ics of the nuclear surface. For this reason we shall study 
the response to the operator dU(r)/drYLM , for L=2 and 
3, where U(r) represents the HF potential calculated in 
the WS cell (cf. Fig. [l| . This operator is specifically sen- 
sitive to the nuclear surface. In the case of an isolated 
nucleus it leads to results similar to those obtained with 
the t l Ylm operator for both quadrupole and octupole 
modes. In the following we shall give a detailed discus- 
sion of the results obtained for the 1364 Sn WS cell. Wc 
shall then present the results associated with the 498 Zr 
case. Finally, we shall consider the evolution of the re- 
sponse in the Sn isotopes with increasing number of neu- 
trons, going from the closed-shell nucleus 132 Sn to the 
drip-line nucleus 1 Sn and then into the inner crust to 
the WS cells 506 Sn and 1364 Sn. 



(a) The case of 1364 Sn 

In Fig. [2] we show the HF and RPA strength functions 
of the operator dU/drYLM in the 2 + and 3~ channels for 
the 1364 Sn WS cell. The bars indicate the solution of the 
discrete RPA equations, while the continuous lines are 
obtained folding the discrete response with a Lorentzian 
function having a Full Width at Half Maximum (FWHM) 
equal to 1 MeV. 

In order to satisfy the Energy Weighted Sum Rule 
(EWSR) within 93% we had to include single-particle 
states up to E cut = 90 MeV. The \ph) basis is composed 
by about 8000 states for the quadrupole response and 
by about 10000 for the octupole response; these values 
are close to the limits of our computational possibilities. 
The fraction of the EWSR obtained by integrating the 
calculated quadrupole response up to a given energy is 



shown in Fig. [3] as a function of energy. We notice that 
the main features of the low-lying part of the response 
can be obtained with a much lower value of E cut . This 
can be seen in Fig. HJa), where the dotted line shows the 
RPA response obtained with E cut = 30 MeV. 

The proton unperturbed response involves deeply 
bound orbitals. The strongest particle-hole transitions 
connect orbitals separated by two oscillator shells, and 
have an energy of approximately 14 MeV, close to the 
value of 2hujp estimated above. There is also a rather 
strong low- lying transition, gg — > ^5/2 > whose energy is 
about 5 MeV. 

The energies and transition strengths associated with 
the strongest (unperturbed) neutron particle-hole transi- 
tions are collected in Tables I (quadrupole response) and 
II (octupole response) . They involve either bound states 
(Eij < E cont — —5.5 MeV) or states in the continuum 
(Eij >= E cont ). The latter are found to be resonant 
states, that is quasibound levels whose wavefunction is 
concentrated in the interior of the nucleus. To show this 
we take a potential which is equal to the mean-field po- 
tential up to 30 fm and is constant and equal to E cont 
from 30 fm to infinity. We then compute the resonant 
states associated to this potential, looking at the wave- 
function phase shift that occurs when a particle at a given 
energy goes past the potential well. The quantum num- 
bers {l,j} label a resonance, if the associated phase shift 
b~ij{E) reaches the value (n + 1)tt/2 (n = 0,1,..) with 
a positive slope at the energy E res . The width of the 
resonance is taken to be equal to the full width at half 
maximum (FWHM) of the derivative of Sij ■ For nar- 
row resonances, the resulting value is almost the same as 
that obtained making use of the expression 

2 

^ = (dS l3 /dE) res - (2) 

Comparing Table U and Table [ED with Table HIl one 
realizes that the quantum numbers of the particle states 
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E (MeV) E (MeV) 

FIG. 2. Quadrupole (a) and octupole (b) strength functions calculated in the 1364 Sn WS cell. Black histograms refer to the 
discrete RPA response, while red and blue histograms refer to the discrete HF response for neutron and protons respectively 
(in units of MeV 2 fm -2 ). The solid curve (in units of MeV fm -2 ) is obtained by a convolution of the discrete RPA strength 
with a Lorentzian curve having a FWHM equal to 1 MeV. The dashed and dash-dotted curve are obtained convoluting the 
neutron and proton discrete HF strength respectively. The dotted line in (a) shows the quadrupole RPA response obtained 
with E cu t = 30 MeV (see text for details). 
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FIG. 3. Fraction of the quadrupole EWSR obtained integrat- 
ing the RPA strength function shown in Fig. HJa) up to a 
given energy. 



coincide with one of those listed in Table Mil (with the 
exception of the state (I, j) = (9, 17/2) in TableHH. The 
wavefunctions of the resonant states in the interior of the 
nucleus are similar to those of bound states belonging to 
the oscillator shells N = 6—10. Their angular momentum 
is the largest available in each shell, except for the shell 
N = 6, where besides the state I = 6, j = 11/2 one 
also finds the states I = 4, j = 7/2 and I = 4, J = 9/2. 
For the highest shells (N = 9, 10) only the state with 
I = j + 1/2 has a resonant character, while its spin-orbit 
partner I = j — 1/2 lies at a too high energy respect to the 
centrifugal barrier. The separation between resonances 
can be estimated from Table Hill to be fc„ « 8 — 9 MeV, 
close to the value estimated above based on the shape 
and the effective mass of the HF potential. 

The value of the spin-orbit splitting for these large val- 
ues of / is close to fvjj n . As a consequence, the intruder 
state (Z, j+1/2) lies rather close to the state (I— 1, j — 1/2), 
and the neutron part of the unperturbed response is char- 
acterized by strong peaks corresponding to transition en- 
ergies of the order of 2hw n (quadrupole response) or of 
the order of lhui n and 2thiu n (octupole response), with 



a spreading caused by the width of the resonances. In 
other words, the strength function is determined to a 
large extent by the shell structure, as in atomic nuclei. 
This point will be further discussed below (cf. Figs. [TT1 
and [P2"|) . 

The RPA 2 + response, shown in Fig. [Ha), displays two 
main peaks, lying close to 3.5 MeV and to 10.0 MeV, with 
a width of about 3.3 and 4.0 MeV respectively. The 3~ 
response (cf. Fig. &b)) shows a strong peak around 3.6 
MeV with a width of about 2 MeV and a broader high- 
energy bump around 18.6 MeV with a width of about 4 
MeV. The neutron and proton transition densities asso- 
ciated with the four RPA solutions carrying the largest 
quadrupole strength are shown in Fig. 01 They are 
peaked at the surface of the nucleus, and their shape 
is close to the derivative of the density, as predicted by 
the semiclassical model for surface collective modes. 

The key role played by resonant transitions can be 
demonstrated by performing a RPA calculation in a WS 
cell with the same number of protons and with approxi- 
mately the same asymptotic neutron density, but with a 
reduced value of Rbox = 25 fm. The resulting strength, 
shown in Fig. [5] is very similar to that obtained in 1364 Sn. 

In order to study in more detail the resonant neutron 
particle-hole transitions in the RPA response, it is useful 
to distinguish their contributions to the strength from 
those involving the other, 'continuum' states [31]. We 
shall limit ourselves to a simple analysis and consider a 
single-particle state in the WS cell to be 'resonant' if its 
quantum numbers coincide with one of those listed 
in Table llnl and if its energy lies in the range E(? s ±.T r ^ s . 
The strength of multipolarity L associated with a given 
phonon \v > can be schematically written 



si = \J2( x ph + (- 1 ) LY p h ) 

ph 

< 3P\\Y L \\]h >< R p \dU/dr\R h > | 2 , (3) 
where X^ h and Y£ h denote the forwardsgoing and back- 
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TABLE I. List of the eight unperturbed neutron particle-hole 
transitions in the 1364 WS cell and associated with the largest 
transition strengths calculated with the operator dU / drY-zu ■ 
In the first four columns we give principal quantum number, 
the orbital angular momentum, the total angular momentum 
and the energies nh,lh,jh and Eh of the hole; in the next 
four columns, the corresponding quantities n p ,l p ,j p and E p 
for the particle. In the last two columns, we give the energy 
Eph — E p — Eh and the transition strength T associated with 
the transition. All energies are in MeV. 



wardsgoing amplitudes associated with each RPA root, 
while R p ,Rh denote the single-particle radial wavefunc- 
tions. In the following, we shall single out the 'reso- 
nant' contributions to the computed strength function, 
restricting the sum in Eq. ([3]) to transitions between two 
bound states (including both protons and neutrons), and 
to 'resonant' neutron transitions, namely those for which 
at least one of the two states has a 'resonant' charac- 
ter. In Fig. [6] we compare the full 2 + and 3~ strengths 
with the correspondent resonant contributions. It is seen 
that the latter reproduce quite well the shape of the full 
response. 

It is also interesting to perform a new RPA calculation 
excluding the transitions between bound states and the 
resonant transitions from the diagonalization. This pro- 
duces a featureless response which accounts for about 5% 
of the EWSR 2+ in the range 0-20 MeV (cf. the dotted 
lines in the insets of Fig. [5]). Performing instead a new 
RPA calculation including in the diagonalization only the 
transitions between bound states and the resonant tran- 
sitions one finds that the resulting quadrupole strength 
displays a very collective peak at E k, 12 MeV with a 
narrow width of the order of 2 MeV, which accounts for 
about 30% of the EWSR (cf. the dashed line in the inset 
of Fig. EKa)). A small peak around 4.5 MeV is also visi- 
ble. This originates from low-energy transitions between 
the well-bound proton states. The octupole strength in- 
stead displays two peaks of comparable strength located 
at about 6 MeV and 20 MeV. The coupling with contin- 
uum states increases the strength of the peaks and lowers 
their energy, producing the full strength: in macroscopic 
terms, this can be associated with a decrease of the sur- 
face tension, taking place when the nucleus is immersed 
in the neutron fluid, as well as with the increase of the 
mass of the fluid which takes part in the motion [32l | . 

From the analysis described above, we conclude that 
the collective response of the system is largely driven by 
the transitions involving the deeply bound protons and 
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TABLE II. The same as in Table [T] for the operator 
dU/drY S M- 
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TABLE III. Total and orbital angular momentum Ij, energies 
E re3 and widths r res of the resonant neutron single particle 
states calculated extrapolating the HF potential of the 1364 Sn 
WS cell up to infinity. The continuum spectrum starts at 
E cont = -3.5 MeV. 



the resonant neutron levels. However, the coupling with 
the continuum shifts the energy of the peaks to lower en- 
ergy and increases considerably the width of the states, 
as well as the strength of the low-energy peaks. We re- 
mark that the low-lying part of the spectrum would be 
modified by neutron pair correlations, not taken into ac- 
count in the present study, although there are no resonant 
single-particle states lying very close to the Fermi energy 
(E F = 5.5 MeV), cf. Table |lri] 
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FIG. 4. Neutron (left panel) and proton (right) transition densities associated with the four strongest RPA transitions calculated 
in the response to the operator dU jdrYiM for the 1364 Sn WS cell, shown in Fig. [2ja). The energies of the transitions are listed 
in the legends of the figures. 
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FIG. 5. The HF (dashed line) and RPA (solid line) quadrupole strength functions calculated in 1364 Sn, are compared to 
those calculated in a WS cell with 50 protons and 654 neutrons and a radius Rbox =25 fm, corresponding to about the same 
asymptotic neutron density (dashed and dash-dotted lines). 
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FIG. 6. (a) The RPA quadrupole strength function calculated for the 1364 Sn WS cell (solid line, cf. Fig[2ja)), is compared 
with the strength obtained including only the amplitudes associated with proton and resonant neutron particle-hole transitions 
(dashed line, cf. text for details). In the inset, the full RPA quadrupole strength function is instead compared with the RPA 
strength calculated using a restricted particle-hole basis limited either to proton states and to bound and resonant neutron 
states (dashed line), or to neutron continuum states (dotted line) (b). The same as (a), for the operator dU / drY^M ■ 
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-7.6 


9 


8 17/2 17.3 


24.9 


94.6 


1 


4 


7/2 


-12.3 


7 


6 11/2 9.0 


21.3 


79.8 


1 


5 9/2 


-1.3 


6 


6 11/2 8.0 


9.3 


55.1 


1 


5 


11/2 


-7.6 


7 


7 15/2 10.1 


17.7 


198.6 


1 


5 9/2 


-1.3 


7 


6 11/2 9.0 


10.3 


73.7 


1 


5 


9/2 


-1.3 


10 


7 13/2 18.6 


19.9 


69.7 


2 


3 5/2 


-0.5 


6 


6 11/2 8.0 


8.5 


35.5 


1 


5 


9/2 


-1.3 


11 


7 13/2 21.4 


22.7 


59.2 


2 


3 5/2 


-0.5 


7 


6 11/2 9.0 


9.5 


42.7 


2 


6 


13/2 


1.5 


9 


8 17/2 17.3 


15.8 


101.4 


2 


6 13/2 


1.5 


6 


7 15/2 9.1 


7.6 


42.0 


2 


6 


13/2 


1.5 


8 


8 17/2 19.4 


17.9 


120.9 


2 


6 13/2 


1.5 


8 


7 15/2 10.1 


8.6 


159.5 


2 


6 


13/2 


1.5 


8 


8 17/2 22.2 


20.7 


68.1 


2 


6 13/2 


1.5 


11 


9 19/2 25.2 


23.7 


75.8 



TABLE IV. List of the eight unperturbed neutron particle- 
hole transitions calculated in the 498 Zr WS cell and associ- 
ated with the largest transition strengths with the operator 
dU jdrYiM- In the first four columns we give the orbital angu- 
lar momentum, the total angular momentum and the energies 
lh, jh and Eh of the hole; in the next four columns, the corre- 
sponding quantities l p ,j p and E p for the particle. In the last 
two columns, we give the energy of the particle-hole jump 
Eph — E p — Eh and its transition strength T. All energies are 
in MeV. 



(b) The case of 498 Zr 

In this section we shall discuss the results of a study 
analogous to that carried out for 1364 Sn, but this time 
for the 498 Zr WS cell. The HF and RPA quadrupole and 
octupolc strength functions are shown in Fig. [7] The 
main differences as compared to 1364 Sn are found in the 
low-energy part of the response, where one notices the 
absence of the low-lying peak in the quadrupole strength 
and the presence of a strong, very low-energy peak in the 
octupole response. Both these features are related to the 
change occurring in the proton single-particle spectrum 
going from Z = 50 to Z = 40. In fact in 1364 Sn one 
finds a gap of about 3 MeV between the last occupied 
proton orbit l<?g/2 and the first unoccupied orbit lgr/2- 
Low-lying quadrupole transitions are then possible be- 
tween the I39/2 orbital and the 1^7/2 or 2 d 5 / 2 orbitals. 
On the other hand, the lowest octupole transitions have 
an energy of about 6 MeV (lftw p ). In 498 Zr the I39/2 or- 
bital is not occupied and low-energy octupole transitions 
are possible (P3/2 — > 39/2 with an energy of 2.7 MeV), 
and only 2 hijp quadrupole transitions are possible. The 
application of these results to the actual case of the in- 
ner crust is model dependent. In fact, as we indicated 
above, detailed calculations of the isotopic composition 
of the crust predict different values for the favoured pro- 
ton number, typical values being Z w 40 or Z w 50. In 
Figs. [TTT e).(f) below we show the strength functions cal- 
culated for 506 Sn (Z = 50). As expected, in this case one 
recovers the low-lying peak in the quadrupole response 
while the low-lying peak in the octupole response disap- 
pears. 

Proton pairing correlations would lead to the partial 
occupation of the lgg/2 orbital, allowing the simultane- 
ous presence of quadrupole and octupole low energy tran- 
sitions. The low-lying part of the spectrum might also 
be affected by neutron pairing correlations. 

The analysis of the neutron transitions leads to the 



TABLE V. The same as in Table [TV] for the operator 
dU/drYsM. 



I 


j 


E res [MeV] 


r res [MeV] 


4 


9/2 


4.7 


2.0 


6 


13/2 


1.6 


5 xlO" 3 


6 


11/2 


8.8 


1.2 


7 


15/2 


10.1 


1.0 


7 


13/2 


21.0 


6.5 


8 


17/2 


19.0 


3.7 


9 


19/2 


30.3 


8.5 



TABLE VI. Total and orbital angular momentum Ij, energies 
Eres and widths L re s of the resonant neutron single parti- 
cle states calculated in the 498 Zr WS cell. The continuum 
spectrum starts at E cont = -0.5 MeV. 



same conclusions presented above in the case of 1364 Sn. 
On the one hand, the transition densities associated with 
the strongest quadrupole phonons (cf. Fig. [8]) have a very 
clear surface character. On the other hand, the strongest 
neutron particle-hole transitions (listed in Tables QV] and 
[Vj are again associated with the resonant single-particle 
states identified by the phase shifts and listed in Table 

Ell 

Furthermore, the analysis of the strength functions 
(cf. Fig. shows that the shape of the response is 
well reproduced including only the amplitudes over pro- 
ton and resonant states. Similarly to 1364 Sn, also in 
498 Zr a RPA calculation including only proton and res- 
onant states produces a narrow collective peak in the 
quadrupole strength, close to 13 MeV (cf. the dashed 
curve in the inset of Fig. [9). Compared to the 1364 Sn, 
however, the coupling to the neutron continuum produces 
a smaller shift of the peak and a smaller increase of the 
width, due to the reduced density of the system. 
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E (MeV) E (MeV) 

FIG. 7. Quadrupole (a) and octupole (b) strength functions calculated in the 498 Zr WS cell. Black histograms refer to the 
discrete RPA response, while red and blue histograms refer to the discrete HF response for neutron and protons respectively 
(in units of MeV 2 fm" 2 ). The solid curve (in units of MeV fm 2 ) is obtained by a convolution of the discrete RPA strength 
with a Lorentzian function having a FWHM equal to 1 MeV. The dashed and dash-dotted curve are obtained convoluting the 
neutron and proton discrete HF strength respectively. 
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FIG. 8. Neutron (left panel) and proton (right) transition densities associated with the four strongest RPA transitions calculated 
in the response to the operator dU IdrYim for the 498 Zr WS cell, shown in Fig. [7J a )- The energies of the transitions are listed 
in the legends of the figures. 
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FIG. 9. (a) The RPA quadrupole strength function calculated for the 498 Zr WS cell (solid line, cf. Fig[7Ja)), is compared 
with the strength obtained including only the amplitudes associated with proton and resonant neutron particle-hole transitions 
(dashed line, cf. text for details). In the inset, the RPA quadrupole strength function is instead compared with the RPA 
strength calculated using a restricted particle-hole basis limited either to proton states and to resonant neutron states (dashed 
line) or to continuum neutron states (dotted line), (b) The same as (a), for the operator dU / drY$M ■ 



(c) Evolution of the response with the neutron number 

Having seen that the resonant single-particle states 
largely determine the 2 + and 3 _ response to the 
dU jdrY^M operator, one may expect that the general 
features of the response to the dU/dr operator should 
evolve in a reasonably continuous fashion going from the 
atomic nucleus into the inner crust. In Fig. [12] we com- 
pare in a schematic way the quadrupole and octupole 
neutron particle-hole transitions calculated 1364 Sn and 
498 Zr with those calculated in the closed-shell, neutron- 
rich nucleus 132 Sn. The horizontal lines indicate the ener- 
gies of the neutron bound levels or of the resonant levels 
(cf. Tables |m] and EJ for 1364 Sn and 498 Zr). The solid 
and dashed lines indicate the strongest quadrupole and 
octupole transitions connecting bound or resonant lev- 
els (cf. Tables U and [ED for 1364 Sn and Tables |IV] and 
[V] for 498 Zr). One can clearly recognize that the shell 
structure plays a similar role in the three cases, leading 
to unperturbed particle-hole transitions characterized by 
an energy 2huj n (for quadrupole transitions) or lhuj n and 
3Huj n (for octupole transitions). 

In Fig. [TT] we compare the the 2 + and 3~ strength 
functions in the inner crust with those of atomic nuclei. 
We keep the same number of protons (Z = 50) and in- 
crease the number of neutrons, going from the neutron- 
rich, closed-shell nucleus 132 Sn (panels (a),(b)) to the 
drip-line, closed-shell nucleus 176 Sn ((c), (d)) and then to 
the WS cell 506 Sn ((e), (f)) and finally to 1364 Sn ((g), (h)). 
The calculation of 506 Sn is very similar to that discussed 
above for 498 Zr, the only difference being the number of 
protons. The proton and neutron mean-field potentials 
are compared in Fig. 1101 

One can see from Fig. [TT] that the main features of 
the giant resonances are the same in the four cases. The 
centroids of the main peaks, associated with 2frw transi- 
tions in the quadrupole response and with ltuo and 3hu 
in the octupole response, are lowered going from 132 Sn 
to 1364 Sn. This is in part due to the shift of the un- 
perturbed response, associated with the larger radii of 
the potentials, particularly in the case of protons, and in 
part is caused by the coupling with the neutron contin- 
uum, as we discussed above. The coupling also increases 
the width of the peaks. A more complete study should 
include the coupling to 2p — 2h configurations. 

In Fig. [UJa),(b) we show also the strength functions 
obtained with the operators rljM and r 3 Y^M in the 
case of 132 Sn. Their shape essentially coincide with those 
calculated making use of the operator dU / drY^M] this is 
not the case for the calculations performed in the inner 
crust, as discussed in Appendix B. 
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FIG. 10. Neutron (a) and proton (b) mean-field potentials calculated for the nuclei 132 Sn and 176 Sn and for the WS cells 506 Sn 
and 1364 Sn. 
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FIG. 11. Evolution of the quadrupole and octupole response going from the neutron-rich nucleus 132 Sn to the drip line and 
into the inner crust. We show the strength functions of the operators dU/drYzM and dU/drY3M calculated in 132 Sn (panels 
(a) and (b)) , 176 Sn (panels (c) and (d)), 506 Sn (panels (e) and (f)) and in 1364 Sn (panels (g) and (h), already shown in Fig. 
[2}. In all the panels, the black histograms refer to the discrete RPA response while the red and blue histograms refer to the 
discrete HF response for neutrons and protons respectively (in units of MeV 2 fm -2 ). The red dashed, blue dash-dotted and 
solid curves (in units of MeV fm~ 2 ) are obtained convoluting the neutron and proton HF responses and the RPA response 
with a Lorentzian function having a FHWM equal to 1 MeV. In panels (a) and (b) we also show the quadrupole and octupole 



responses associated with the operators r 
and multiplied by 120). 



'Yim (in units of fm and multiplied by 4) and r Y aM in Sn (in units of fm 
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FIG. 12. The strongest neutron particle-hole transitions in 132 Sn (a), in 498 Zr (b) and in 1364 Sn (c), see text for details. 
Single-particle levels of even and odd parity are drawn by solid and dashed lines. Solid (red) and dashed (blue) arrows refer to 
2 + and to 3~ transitions. The Fermi energy is represented by the thick dashed line. 
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IV. CONCLUSIONS 

As Negele and Vautherin noticed in their seminal pa- 
per, 'the degree to which the nuclei in the free neutron 
regime resemble ordinary nuclei' is striking, and 'this sim- 
ilarity is also manifested in the behaviour of the single- 
particle energies' 251 ] . In this work we have shown that 
the persistence of the shell structure influences the lin- 
ear response of the Wigner-Seitz cell in the inner crust 
to a large extent. We have found it useful to study the 
response to the operator dU/drYLM, which acts as a fil- 
ter, enabling one to focus on the effects associated with 
the nuclear surface. Due to the existence of resonant 
states with rather narrow widths, the main features of 
quadrupole and octupole giant resonances are similar to 
those of ordinary atomic nuclei. However, the interac- 
tion of the bound nucleons with the neutron sea, has an 
important effect, leading to a decrease of the energy of 
the peaks, and to an increase of their width. 

The existence of vibrational modes associated with the 
nuclear surface can have important consequences on the 
neutron superfluidity of the system in the So channel. 
While medium effects tend to suppress pairing correla- 
tions in neutron matter, being dominated by spin fluc- 
tuations, the surface fluctuations studied in this paper 
lead to an attractive contribution. Even if the global ef- 
fects may not be large, due to small volume occupied by 
the nucleus, they can affect the spatial dependence of the 
pairing gap, and have important consequences on vortex 
pinning [If @, In this respect, it would be also im- 
portant to examine the coupling of neutrons with lattice 
vibrations 1351. [3611. 
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VI. APPENDIX A 

In this Appendix we shall compare the results of nu- 
merical calculations of the linear response in a WS cell 
without protons with analytical results which are avail- 
able in uniform matter in the case of Skyrme interactions 

& 

We shall consider the excitation of either density or 
spin modes, computing the response function to the ex- 
ternal fields 



uf-r 



6 («) 



(4) 



where a = 0, 1 distinguishes the two possible channels 
5 = and 5=1: 0° = 1, O 1 = er 2 according to whether 
density or spin modes are considered. 

In neutron matter, the RPA strength function can be 
written as 



S {a) (q,E) = --ImIL (a \q,E) 



-- J d 6 rd 6 r'V ext {?)V ex t{7)Imn{r, ?, E), (5) 

where II(r, r 1 , E) is the polarization function [34| . 

The RPA strength function calculated in neutron mat- 
ter with a Skyrme effective interaction reads 



n^(v,k) = 2n [i-iu 1 (Q) n 

k 2 - 



2wt ] k 2 F 



l-(m k k 3 f /3-K 2 )W^° 



fin 



(w^ a) k 2 f y (n2-n n 



4k 2 v 2 n 2 



2W 2 (a) fc 2 f (2fc 2 n -n 2 



(6) 



where H2 e II4 are generalized Lindhard functions given 
and we have introduced the adimensional vari- 



ables k = jjL- and v — ff*^ , expressed in terms of the 
Fermi momentum kF and of the effective mass mu- The 
coefficients w\ can be expressed in terms of the usual 
coefficients defining the Skyrme interactions. In neutron 
matter one finds 



wl = 
w\ = 



(7 + l)(7 + 2) 



1 3 
t Q {l - x ) + -ii9 2 (l - xi) - -t 2 q 2 {l + x 2 ) 

-t (l ~ x ) - \hq 2 (l - si) - ^t 2 q 2 (l + x 2 ) 
*3 



(. 



P 7 (l-*3) 



|ti(l-ari) + -t 2 (l + a!2) 
--t 1 (l-x 1 ) + -t 2 {l + x 2 ) 



(7) 



The unperturbed response is obtained substituting the 
Lindhard function ILj in place of Tl(q, E) . 

We shall compare these analytic results with the nu- 
merical calculations performed putting 498 neutrons in a 
WS cell of radius Rws = 42.2 fm without protons. In 
the latter case the response function reads 

+00 

S { wl(9,E) = V \<n I e^eW | >\ 2 L(E, E n ), (8) 



n=0 



where |0) is the ground-state, \n) are the (mean-field or 
RPA) excited states of the system, and the computed 
discrete response has been convoluted with a Lorentzian 
function L(E, E n ). 
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Exploiting the spherical symmetry of the system, one 
can introduce the multipole decomposition of the exter- 
nal operator 



+oc +L 

-tvrV V t L j L (qr)Y£ M (n q )Y LM (n r ), (9) 



L=0M=-L 



where jl (qr) is a spherical Bessel function, leading to an 
analogous decomposition of the response in the WS cell: 

+00 

The RPA excited states are given by 
I"-} = \vjm) = 

E Ph (xZk(J)Wh)- 1 ) + Y; h (J)(-y +M \( P )- 1 h}) , 

( 10 ) 

where the X and Y denote the forward and backward 
amplitudes calculated in RPA, and one obtains for den- 
sity modes (S = 0, J = L) 

sWs.i.ti, E ) = Mil + 1) E, =0 E„ h (x; h (L)+ 



while the RPA response function in the S = channel is 
given by 



fEWSRfr 



RPA 



V 



1„2 



unif 



Kg 
' 2m 



(14) 




FIG. 13. Response function per unit volume as a function 
of the energy in the 5 = channel, computed in a WS cell 
of 42.2 fm radius containing 498 neutrons. The contributions 
from several multipolarities are shown. The linear momentum 
transferred by the external field is q = 0.51 fm . 



(-) L Y1(L)) (j p \\i L Y L \\j h )(R p \\j L (qr)\\R h ) 



L{E,E V ) 
(11) 

The contributions to the S = response in the cell as- 
sociated with several multipolarities are shown in Fig. 
I13[ for the momentum transfer q = 0.51 fm -1 . It 
is seen that the contributions from the different mul- 
tipolarities increase with L reaching the maximum for 
L = 15, and then rapidly decrease, becoming negligi- 
ble for L larger than about 20. This is in keeping with 
the simple estimate of the maximum angular momentum 
which can be transferred by the external held, given by 
L m ax ~ qRws ~ 21. 

The total response S ws (q, E) to the external fields (j4]) 
is compared in Fig. Q3]with the analytic result, in uni- 
form neutron matter at the same density (p ~ O.Olpo) 
and for the same value of q =0.51 fm -1 , both in the 
unperturbed and in the RPA case. We find reasonable 
agreement with the calculations in the WS cell. The 
strength functions computed in the WS cell display a tail, 
which is not present in the uniform system, where the ex- 
citation energy of the system has a maximum value, im- 
posed by energy and momentum conservation and equal 
to 



h 2 



2mi 



{q A + 2qk f ), 



(12) 



which, in the present case, corresponds to about 13.5 
MeV, taking into account that the neutron effective mass 
is very close to the bare mass (m^ = 0.99 m) at the 
present density. The analytic Energy Weighted Sum Rule 
(EWSR) associated with the unperturbed response func- 
tion S(q, E) is given by 



EWSR H f\ 



2„2 



V 



unif 



} , 

Irak 



(13) 




FIG. 14. Response functions per unit volume as a function of 
the energy in the 5* = and 5 = 1 channel, computed in a 
WS cell of 42.2 fm radius cell containing 498 neutrons (dashed 
lines) and in uniform neutron matter (solid lines) for the linear 
momentum q = 0.51 fm - transferred from the external field. 
Both systems correspond to a Fermi momentum fcj? = 0.39 
fm" 1 . 

The Energy Weighted Sum Rule (EWSR) in the WS 
cell is numerically calculated as the sum of the EWSR 
associated to different multipolarities J: 

EWSR (e i?? ) = 4ir^2EWSR{j L (qr)). (15) 

L 

The contribution calculated for each L can be compared 
to the value obtained evaluating the classical EWSR on 
the HF ground state: 



EWSR c i as 



i(jL(qr)) 



£;^N < 0| ({iMqr)) 2 + L(L + 1) (lj L (qr)) 2 ) |0 >, 

(16) 
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The cumulative EWSR is shown in Fig. [15] for the RPA 
calculation. The total value and turns out to be only a 
few per cent larger than the analytic value. 
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FIG. 15. EWSR per unit volume calculated in the 5 = 
channel for as WS cell of 42.2 fm radius containing 508 neu- 
trons, which corresponds to a Fermi energy Rf = 0.39 fm -1 . 
The linear momentum transferred from the external field q — 
0.51 fin -1 . The squares show the cumulative EWSR obtained 
from the RPA strength function as a function of the orbital 
angular momentum L, while the dots are obtained using Eq. 
(|16|) . The arrow indicates the analytic value for the uniform 
system (cf. Eq. (O). 

We conclude that the numerically computed response 
in the WS cell without protons is quite close to the re- 
sponse calculated analytically in uniform neutron matter 
at the same density. This gives us confidence in our nu- 
merical approach, which is applied in the main text to 
the case of a WS cell with a nucleus in its center. 



show the features discussed in refs. [21] and (22[, dis- 
playing a strong bump at low energy ('named supergiant 
resonance' in refs. [Il|,[22j]). The presence of the nucleus 
does not change the energy dependence of the strength, 
except for the fine details of the fragmentation of the 
discrete peaks produced by Landau damping. The main 
effect of the nucleus is the decrease of the absolute value 
of the EWSR: the EWSR for 1314 X is about 20% smaller 
than for 1364 Sn, due to the fact that in the latter about 
150 neutrons are bound in the nucleus, and couple to the 
r 2 operator much less efficiently than the unbound neu- 
trons. The neutron transition densities associated with 
the vibrations bulding the main peak are shown in Fig. 
El for the 1314 X and for the 1364 Sn cell. In both cases 
they show large values in the middle of the cell. The nu- 
cleus in the 1364 Sn case only induces a modest increase 
of the transition densities for r ~ 10 fm for two of the 
phonons. These results are obviously in keeping with the 
fact that the r 2 operator weights heavily the region far 
from the nucleus, at variance from the dU / drY^M oper- 
ator considered in the main text. 



VII. APPENDIX B 

In ref . [U [22j , the strength associated with the op- 
erator r L was calculated in the inner crust in RPA and 
in QRPA (see also [20(). Usually the operator r L is used 
as a convenient approximation to the operator e lq ' r in 
the limit of long wavelengths. In the inner crust one is 
dealing with boxes with of the order of 20 fm and with 
Fermi energies of the order of several MeV, so that the 
produce kpR is not small and the choice of this operator 
does not seem to be well motivated. 

Nevertheless, we have calculated the strength function 
associated with r 2 Y2M for the cells 1364 Sn and 498 Zr in 
order to compare with the results obtained in [22|. The 
strength functions are shown in Fig. 1161 where we also 
show calculations for two cells without protons and with a 
similar number of nucleons, indicated as 1314 X and 508 X. 
In the present case, a cut off E cut — 30 MeV is sufficient 
to exhaust the EWSR. 

Our results coincide with those calculated by the Orsay 
group for the same cells [l?]], and the strength function 
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FIG. 16. (top) HF and RPA quadrupole strength functions calculated in the cells 1364 Sn (a) and 1314 X (b), in units of MeV" 1 
fm 4 . (bottom) The same but for the cell 498 Zr (c) and 500 X (d). The transition strengths of the discrete states calculated in 
RPA are also shown in histogram form (in units of fm 4 ). 




FIG. 17. Neutron transition densities associated with the strongest phonons calculated in 1314 X (a) and in 1364 Sn (b). The 
energy of the phonons is indicated. 
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